Formalism for obtaining nuclear momentum distributions by the Deep Inelastic 

Neutron Scattering technique 
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We present a new formalism to obtain momentum distributions in condensed matter from Neutron 
Compton Profiles measured by the Deep Inelastic Neutron Scattering technique. The formalism 
describes exactly the Neutron Compton Profiles as an integral in the momentum variable y. As 
a result we obtain a Volterra equation of the first kind that relates the experimentally measured 
magnitude with the momentum distributions of the nuclei in the sample. The integration kernel 
is related with the incident neutron spectrum, the total cross section of the filter analyzer and the 
detectors efficiency function. A comparison of the present formalism with the customarily employed 
approximation based on a convolution of the momentum distribution with a resolution function is 
presented. We describe the inaccuracies that the use of this approximation produces, and propose 
a new data treatment procedure based on the present formalism. 



I. INTRODUCTION 

During the last decades there has been a growing in- 
terest in the experimental determination of momentum 
distributions of different particles in a wide variety of 
systems, being of principal importance the experiments 
based on deep-inelastic scattering techniques. Thus, the 
Compton scattering of x rays, gamma rays, and electrons 
has been for a long time used to obtain the electron mo- 
mentum distribution in condensed matter, the momen- 
tum distributions of nucleons in nuclei have been deter- 
mined by inelastic scattering of high-energy protons and 
electrons, and the deep inelastic scattering of electrons, 
muons and neutrinos was essential for establishing the 
quark model of the nucleonsi^. 

The experimental determination of atomic momentum 
distribution in condensed matter making use of differ- 
ent neutron techniques has experienced a considerable 
evolution up to the present times. An historical review 
of these neutron experimental techniques has been pre- 
sented in Ref. 0- Notably, the Deep Inelastic Neutron 
Scattering (DINS) technique, proposed by Hohenberg 
and Platzmann 4 - nearly 40 years ago for the study of the 
Bose condensation in superfluid 4 He, attained a consid- 
erable experimental development. The technique proved 
to be remarkably adequate to explore the momentum 
distributions of the nuclei in condensed mattes£*&, and 
its development continues nowadays^. Being it a tech- 
nique based on large momentum and energy transfers, 
the interpretation of the recorded spectra is based on 
the validity of the Impulse Approximation (IA), where 
it is assumed that the target particle recoils freely after 
the collision with the incident particle^, assuming that 
the binding energy of the target nuclei is negligible com- 
pared with the energy transferred by the neutron to the 
target nuclei 4 *^^. The justification of the validit y of the 
IA in DINS experiments can be found in Refs. [ll] and 
[Ta, while the deviations either due to final-state effects 



(FSE) 
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, or to initial-state effects (ISE)A£ were 



extensively treated. The IA allows a direct relation of the 
scattering law with the momentum distribution function, 
making use of the variable y, which is the projection of 
the initial impulse of the target nuclei over the direction 
of the momentum transfer vector q, through a process 
known as y - scaling 1 ^ . 

On the other hand, the reactor based techniques for 
the measurement of atomic momentum distributions (de- 
veloped before DINS) have employed triple-axis crystal 
spectrometers, magnetically pulsed beams, or rotating 
crystal and chopper time-of-flight spectrometers^. Since 
these techniques employ neutrons with incident energies 
restricted to the thermal range, they do not satisfy the 
IA and therefore sophisticated corrections by FSE be- 
come essential before the momentum distribution can be 
obtained. The advent of high epithermal neutron fluxes 
available from intense pulsed spallation sources, placed 
the DINS technique as the most suitable one for the study 
of momentum distributions of nuclei constituting the con- 
densed matter. 

A common practice in the data treatment of the DINS 
experiments is to employ a convolution approximation 
(CA)£ii2ii£, wnere it is assumed that the observed spec- 
trum in the time-of-flight scale can be expressed as a 
convolution in the variable y of the momentum distribu- 
tion, with a mass-dependent ad hoc resolution function 
devised to take into account the experimental uncertain- 
ties along with the resonant filter energy width. However, 
the accuracy of this procedure has been questioned by us 
in previous papers22*2i. The criticisms rested on the fact 
that the expression of the CA is not directly deducible 
from the exact one. When we tested it against the exact 
expression we observed inaccurate results for the peak po- 
sitions, areas and shapes of the observed profiles, and the 
effects were more noticeable when treating light nuclei. 
As a consequence, a noteworthy result of DINS reported 
on H2O/D2O mixtures^, that reported to have found an 
anomalous behavior in the neutron cross sections, was 
critically analyzed from the point of view of the accuracy 



of the CA. A careful analysis lead us to conclud o 20 i 21 that 
the inaccuracies committed in the CA produce anomalies 
in the same sense as reported. Later measurements per- 
formed by transmission on H2O/D2O mixtures 2 ^, showed 
no traces of anomalous neutron cross sections. Further- 
more, on the theoretical side, recent publications cast 
doubts on the existence of anomalous phenomena 24 - 25 . 
The controversial situation increases the need to revise 
the procedures on which the CA is based. 

Besides the mentioned limitations, there is another im- 
portant unsatisfactory feature of the CA, that will be 
treated in the present paper: in the CA only a single final 
energy corresponding to the peak of the resonant filter is 
considered when the transformation from time-of-flight 
to the y scale is performed. However, as was remarked 
in a recent analysis^, the final energies of the scattered 
neutrons are not restricted to a narrow resonant-filter line 
width, but to a broader energy distribution determined 
by the kinematics, the dynamics and the full total cross 
section of the filter. The contributions of detected neu- 
trons with final energies outside the main resonance peak 
cannot be neglected and in some cases become dominant. 

The motivation of the this paper is to expose the de- 
tails of the exact treatment that must be performed on 
the DINS experimental data, in view of the inconve- 
niences that stem from the mentioned weakness of the 
CA. We reformulate the basic equations that describe the 
neutron Compton profiles starting from the exact expres- 
sions as was previously treated 2 ^. As a result, explicit in- 
tegrals in the variable y are obtained. The passage from 
the time-of-flight variable t to y is presented in its ex- 
act form, and the final expression results in a Volterra 
equation of the first kind, which can be solved for the 
momentum distribution J(y). The results of the present 
formalism, which is essentially exact are compared with 
those produced by the CA. We stress on the imperfec- 
tions that it generates and how those problems are solved 



by employing the present formalism. The expressions we 
present are amenable for the use of the experimental- 
ists, and we recommend them to replace the currently 
employed formalism in the data analysis of DINS exper- 
iments. 



II. THEORETICAL BACKGROUND 

A. Preliminary considerations 

In this section we will present a brief summary of the 
basic equations that will be employed throughout this 
paper for describing the observed spectra in DINS experi- 
ments. For a more complete theoretical development, the 
reader is referred to Ref. |2{J. We will restrict our anal- 
ysis to a standard inverse-geometry DINS experimental 
setupi&, in which incident neutrons with energy Eq (char- 
acterized by the spectrum $(i?o)) travel along a distance 
Lq from the pulsed source to the sample. Scattered neu- 
trons at an angle 6 and a final energy E travel along a 
distance L\ up to the detector position. A movable filter 
with a neutron absorption resonance in the eV-energy 
region is placed in the scattered neutron path, so con- 
secutive 'filter in' and 'filter out' measurements are per- 
formed. The spectra are recorded as a function of the 
total time of flight t. Throughout this paper, we will ex- 
plicitly omit the description of experimental uncertainties 
due to geometry, time of flight or multiple scattering. For 
an account on these effects the reader is referred to Refs. 
GJandlll 

The measured magnitude (known as Neutron Compton 
Profile (NCP)) is the difference count rate ('filter out' 
minus 'filter in') as a function of t, which for a point-like 
sample can be expressed as2& 
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where a(Eo,E,9) is the sample double-differential cross 
section, e(E) the detector efficiency and AO the solid an- 
gle subtended by the detector. The term between brack- 
ets is the absorption probability of the resonant filter, 
characterized by a number density n, a thickness T and 
a total cross section <jf(E). As indicated, the integral in 
Eq. J3> must be calculated at a constant time-of-flight t 
given by the kinematic condition 



t = 



Lq 
jE~n 



s/E 
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where m is the neutron mass. The lower limit of integra- 



tion -Eoinf is determined from the kinematic condition 
that in the second flight path the neutron has infinite 
velocity. Its value is 



E, 
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As usual, let 

huj = E - E (4) 

be the energy transferred by the neutron to the sample, 



and 



hq = J 2m (e q +E- 2^ E Ecos9 



(5) 



the modulus of the transferred impulse. The DINS tech- 
nique was developed on the theoretical basis of the IA, in 
which the scattering law (for a monatomic sample) can 
be writter 
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where e q is the unit vector along the direction of q. The 
variable y is the projection of the impulse p of the nucleus 
of mass M on e q , and can be expressed as 



V = fq ( 



2M 



(7) 



J(e q , y)dy is defined as the probability that a nucleus has 
a momentum component along the direction of e q with 
values between y and y + dy. The distribution J (e q ,y) 
must be symmetric around y — if a moment in a given 
direction of the space is equally probable than in its oppo- 
site. For an isotropic sample J (e q , y) does not depend on 
the direction given by e q , and then J(e q ,y)dy — J(y)dy 
is the probability that a nucleus of the sample has com- 
ponent of momentum between y and y + dy along any 
direction in the space. So, the dynamic structure factor 
is reduced to 



S(q,w) 



M 

hq 



J(y) 



(8) 



Then, for an isotropic sample composed of different nu- 
clei, being Nm the number of nuclei of mass M and 
bound scattering length bM, the double differential cross 
section is 



a(E , E, 9) = -L^ILj2N M b 2 M MJ M (y M ) 



(9) 



where the sum is extended over all the masses corre- 
sponding to the different species (in non equivalent posi- 
tions) present in the sample. 



B. Convolution approximation 

The usually employed convolution approximation is ex- 
pressed by«i& 

Cconv{t) = ^2 ^Af(*)JAf (j/Af) ® RmO/m), (10) 



M 



with 



&f(t) = N M b M M 
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(11) 



where E\ is fixed, defined by the energy of the main ab- 
sorption peak in the filter total cross section. The sum 
has the same meaning as in Eq. ©, and a resolution 
function R(%im) is introduced as a way to contain the ge- 
ometric uncertainties as well as the filter resonance width. 
It is worth to emphasize that that in Eq. i|l(J|) . the re- 
lation between the variable yM and t is made through 
Eq. Q , where q and u are calculated also with the same 
Ei fixed, and Eq compatible with the time-of-flight rela- 
tion J2}. Also, it must be noticed that the term between 
square brackets in Eq. Ijlll) is independent of the sample 
characteristics. 

It is important to notice that the resolution function 
in the CA framework is deducible from Eq. i|l(J|l and Eq. 
considering a sample represented by an ideal gas in 
the limit of T — ► OK. In this case Jm(Vm) — > S(vm) an< 
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When geometric uncertainties are considered, they also 
contribute to the width of Rm , but as we have already 
mentioned, this case will not be analyzed in the present 
paper. It is worth to comment that in common practice 
Eq. I |12|l is not employed, but instead a fitted Lorentzian 
(for Gold filters) or Gaussian function (for Uranium) in 
the variable y are employed^. 



III. FORMULATION OF THE EXACT 
EXPRESSION 

As it was commented, Eq. (fTUf can not be deduced 
from the exact expression |QJ and therefore its applica- 
tion can not be fully justified. On the other hand an 
expression where (like in l|10|l ) the y variable appears ex- 
plicitly, is desirable since it is the natural variable of the 
impulse distribution J(y). In this section we will deduce 
such exact expression that keeps y as the integration vari- 
able. To this end, we will change the integration variable 
from Eq to y in Eq. ifTJl. In the first step we rewrite Eq. 
in terms of E, Eq and 9 as 
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M [Eq-E [Eo + E-2y/E Ecos9 



2m (Eq + E - 2^Eq~Ecos8) 



(13) 

At a given time of flight t, E and E are linked through 
Eq. J2J. Therefore, at a scattering angle 8, yM is a 
function of Eq and t. On the other hand it must be 
noticed that Eq can be a multi-valuated function of j/m- 
There is a substantial difference between this definition of 
yM and that employed in the CA (Sect. IIIBjI . called yM- 
In the present definition j/m is a function of Eq for each 
considered time of flight, while in the former definition 
yM takes one value for each t. 

To proceed with the variable change, we first study the 
limits of integration in the variable y. When Eq — ■» £t)mf 



(defined in Eq.©), y — > — oo regardless the mass of the 
scattering nuclei. When E — > +00, it can be readily 
shown that 



In Sect. IIV Al we will present a detailed study of the 
behavior of tjm- 



lim y M 

Eq— >+oo 
t,8— const. 



if M > m; 
if M < m; 



m% cos if M = m. 
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The sought variable change in Eq. QJ involves the 
Jacobian (calculated in Appendix lA"|l 
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To perform the variable change from Eq to y in the in- 
tegral of Eq. QJ, it is necessary to know the intervals 
where the function y(E ) is monotonous. Therefore we 
must obtain the Eq values (and its corresponding y val- 
ues) where 



dy_ 

dE 
Defining the variables 



= 0. 



(16) 



(17) 



(where the positive value of the square root is taken in 
the definition of x), Eq. I|ltj|) can be reduced to a fifth- 
order polynomial in x, with real coefficients. The sought 
values of Eq are obtained from the roots of 
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m - 3M + (m + M)^ cos + 2to cos 2 (19) 

(M - to) cos - (3M + to)£ - 2ml cos 2 
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Therefore, Eq. (|18|l has at most five real roots and 
at least one, and only the positive values have physical 
sense. Let Nr be the number of such roots (0 < Nr < 5), 
whose values will be labeled Xj, and correspondingly Eq 
and yu,j the roots in the Eq- and y-scales respectively. 
Now we can proceed to change the integration variable 
from Eq to y in Eq. ©. To this end, the integration 
range has to be partitioned in Nr+1 intervals (Eq~ ,Eq), 
where the function yAf(^o) is monotonous, which we will 
identify as 'branch j' (where E® — E oin f given in Eq.© 
and E^ R+1 = +00). The result is 
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where the superscript j in dEo/dynj indicates the branch where it must be evaluated, yu,a — ~oo and yM,N a +i is 
defined in Eq. g3J. Thus, the NCP results 
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is defined in the interval (jjM,j-i,yM,j), and is zero out- 
side such range. 

Finally, for a sample containing Nm atoms of mass M 
in non equivalent positions, the NCP can be expressed 
as 



=(*) = X> M (t) 



(23) 
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where 



c M (t) 



+ 00 



dVMfM{VM,t)J M {VM) (24) 



and fhiiyM, t) is the sum on the different branches (1 < 
j<N R + l)oif M (y M ,t). 

Eqs. I|21fl and (|22|l are the central expressions of the 
present paper, which are exact expressions, since they 
were derived directly from Eq. Q) without any approx- 
imation. Eq. H21fl is a Volterra equation of the first 
kind. To obtain J(y) from it, it is possible to resort to 
numerical methods. Alternatively, if the function that 
describes J(y) is known beforehand (or a plausible form 
is assumed for it), it will be possible to obtain it from a 
fitting procedure. In this last case it will not be necessary 
to change the integration variable to y. In the case where 
the sample and/or detector dimensions are not negligible, 
Eqs. Ij20l24|l must be integrated also over the geometric 
dimensions. 

As it was already mentioned, the definition of j/m em- 
ployed in Eqs. (|20I24() differs from that employed in the 
CA. It is worth to notice that while in the CA the resolu- 
tion function is the same for every t (centered in yiait)), 
in the present formalism, /m(j/M>£) depends explicitly 
on t. It must be remarked that if the CA framework 
were exact, Eqs. (|10|) and (|21|l would produce the same 
NCP, for an arbitrary J(y). Thus, from a direct compar- 
ison between Eqs. (|24() and (|10fl the following relation 
between the exact expression of fuiyMji) and the reso- 
lution function would be verified 



$(t)R M (v M (t)-y) = f M (v,t) 



(25) 



where £m(£) is a defined in Eq. i|ll|) . However, the usu- 
ally employed Rm(]Jm) does not verify expression (12511 . 
as will be examined in the next section. 





r\ H 




■ 


1.5- 






- 


1.0- 






- 


0.5- 


\ / \ 


D 


- 








■ 


■ 


1 ■ 1 


1 





t(ns) 



FIG. 1: (Color online) Neutron Compton profiles of hydrogen 
and deuterium at a scattering angle of 9 = 70° calculated 
considering as flight paths La = 1105.5 cm and L\ — 69 cm, 
assuming a gold foil analyzer of nT = 4x 10 -5 barn -1 and 
effective temperatures of 115.22 and 80.11 meV respectively. 



described by $(E ) = E^"' 9 , and flight lengths L = 
1105.5 cm and L\ = 69 cm (so £ = 16.022) will be chosen 
to match the layout of the DINS facility at the Ruther- 
ford Appleton Laboratory^. The filter will be represented 
with a gold foil of nT = 4 x 10 -5 barn -1 and its total 
cross section will be described from the data compiled in 
Ref. i3 

We will divide our discussion in three different aspects 
of our formalism. The first will be devoted to a close 
examination of the behavior of the variable y. We will 
establish the relationship between y in the exact formal- 
ism and in the CA framework. In the second part we will 
analyze the behavior of the kernel fmiy^t) and we will 
compare it with the resolution function employed in the 
CA. In the third part we will analyze the consequence of 
our analysis in a specific case. The starting result over 
which we will base the subsequent discussions are the 
NCPs of hydrogen, deuterium (Sects. IIV Al and IIVB|) 
and 3 He (Sect. IIVC|) as calculated from Eq. (JTJ. 

The study of H and D will be based on ideal gases of 
effective temperatures of 115.22 and 80.11 meV respec- 
tively, as described in Ref. |20|. Their NCPs are shown 
in Fig. [1] 



IV. DISCUSSION 

In this Section we will analyze different aspects of the 
formalism we are presenting, with especial emphasis on 
the differences between the exact results and those ob- 
tained in the CA framework. We will focus our attention 
on scattering on H, D and 3 He, since as we will show, 
they are the most sensitive cases in normal experimen- 
tal conditions. We will assume throughout detectors of 
ideal efficiency e(E) — 1. An incident neutron spectrum 



A. Study of the variable yu 

We will begin our discussion by examining the behav- 
ior of the yu variable as a function of the time of flight. 
In Fig. |21 we show its behavior according to Eq. 1131) as 
a function of the incident neutron energy Eq at different 
times of flight, for the case of hydrogen (M/m=0. 99862) 
at 8 =70 ° . The behavior at the limits commented in 
Ea. (|14|) are observed. In the same figure we show the 
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FIG. 2: (Color online) The impulse y as a function of the 
incident energy Eq for A=0. 99862. The exact calculations 
according to Eq. l|lMfl at the indicated times of flight are 
compared with the calculation performed by assuming the 
final energy fixed at 4.906 eV, which is the variable y employed 
in the CA framework 



behavior of yu according to the conditions normally em- 
ployed in the CA framework (Sect. IIIBjI . The differences 
between both ways to define y is evident: while in the ex- 
act formalism there is a whole range of possible values of 
y at each time of flight (which is a function of Eq), in 
the CA there is only one value of y defined at each t. 
We observe the curves y(Eo) in the exact formalism in- 
tersects the one of the CA in two points. However at 
each time of flight there is only one value of y defined in 
the CA framework. Thus, only one of those two points 
over the exact y-curve has a final energy corresponding 
to that of the filter (4.906 eV), while the other has the 
same value of y and -Bo but a different E. At 85.22 /xsec, 
both points coincide so the curves are tangent. At times 
of flight greater than 85.22 /zsec, this property is found 
in the intersection of the leftmost point between both 
curves, while at times lower than 85.22 /xsec we find it in 
the rightmost intersection point. 

The behaviors observed in Fig. |2 are clarified in 
Fig. Ola), where we mapped on the (Eg , E^-plane the 
constant-y and constant-i curves, for the case of hydro- 
gen. In Appendix iBl we show the details for the calcula- 
tion of the constant y curves. We also show the constant- 
E line corresponding to main peak absorption energy of 
the gold analyzer (E=4.906 eV). The intersection of a 
constant-t curve with the 4.906 eV-line defines an Eq- 
value with which the y-value in the CA framework is 
calculated. It can be easily shown that in this case for 
each constant-^ curve, there is one (and only one) tangent 
constant-y curve. The geometric loci of the tangency 
points in the (Eq, E)--p\a,ne result to be straight lines 
(which will be called limit lines (LL) hereafter) whose 
slopes can be obtained from the real positive roots of 




FIG. 3: (Color online) (a) The full thick lines indicate the 
constant-?/ curves for hydrogen at 9 — 70° and £ — 16.022. y- 
values are indicated in A -1 . In dotted lines, the curves of t = 
constant (the indicated values are in /tsec). The horizontal full 
line corresponds to E=4.906 eV. (b) Detail of the upper frame 
showing a case where a constant-i and a constant-y curve are 
tangent. The particular case of a tangency at i?=4.906 eV 
was selected (see text for details). 



Eq. iJTSJ. Bearing in mind Eq. i|17|) those slopes are the 
squares of such roots. In this particular case of hydro- 

dy 



gen 



dEo 



= has only one real root, which defines 



the line E — 0.021^0 as the sought geometric locus of 
the tangency points, which is represented in Fig. Bfb). 
In the representation of Fig. [21 this corresponds to the 
geometric locus of the maxima of the y(Eo) curves. It is 
worth to point out that the consequence of having one 
real root in Eq. (|18|) is that there are two branches, so 
N R = 1 in Eq. J2TJ. In Fi g- Mj°) w e show a detail of the 
upper frame. The line E = 0.021i?o intersects the con- 
stant E = 4.906ey at a point corresponding to t =74.91 
/xsec and y = 9.73 A -1 . In consequence, in Fig. [5] the 
intersection between y and the curve of 74.91 /xsec oc- 
curs at its maximum. The importance of the LL resides 
in that it indicates the y values where /m(z/,£) presents 
singularities, and will be commented in the next section. 
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FIG. 4: (Color online) The open symbols indicate value of 
the slopes of the straight lines corresponding to the geometric 
loci of the tangency points between constant-^ and constant- 
y curves in the (E , £t))-plane as a function of the scattering 
angle, for H, D and O, with a flight-paths ratio £ = 16.022. 
Full symbols indicate the ratio E/Eq in the cold-gas limit for 
H, D and O. 



The slope of the LLs is a function of the scattering 
angle and I. This is shown in Fig. 0J for I = 16.022. In 
the case of hydrogen a LL can be defined in the whole 
angular range from ° to 90 ° existing only one of such 
LLs as already mentioned. On the other hand, for deu- 
terium there exist two real positive solutions of Eq. Ijl8|l 
below 39.7 ° and so two LLs, and for oxygen we observe 
the same situation at a limit angle of 24.7 " . 

It is also interesting to observe in Fig. the represen- 
tation in the (Eq, _E)-plane for deuterium at 9 =35 ° . In 
this case there exist two LLs (in accordance with Fig. 0J, 
which are shown in the graph. 



B. Kernel f M (VM,t) 

Having established the behavior of the variable yu we 
can now analyze Eq. (|24l) . We are especially interested 
in examining the kernel /m(?/m,£)- In Fig. 0we show 
its behavior in the case of scattering in hydrogen at 70 ° 
for different times of flight of interest in the Compton 
profile, as can be checked in Fig. ^ It is very impor- 
tant to notice the asymptotic behavior of /m (yu , t) at 
the values of y^ defined by the LL, that is understood 

dEg 
dy 

(|22|l . In the same figure we show the corresponding ker- 
nel in the CA framework expressed in the left hand side 
of Eq. Ij25[) . where the resolution function was calculated 
according to Ea Q12[l . The maxima of both distributions 
are observed to be located at ijm (t), i.e. the y variable as 



from the fact that the Jacobian 



oo in Eq. 



i.O 



FIG. 5: (Color online) The full thick lines indicate the 
constant-?/ curves for deuterium at 9 = 35° and £ — 16.022. 
y- values are indicated in A -1 . In dotted lines, the curves 
of t = constant (the indicated values are in //sec). The line 
corresponding to y =0 is shown in dash-dotted line. 



defined in the CA framework. However, the asymptotes 
present in the exact formulation mark important differ- 
ences between both formulations, which are clearly man- 
ifested at 74.91 /isec, as showed in Fig. |HJb), where the 
asymptote is exactly in jJm (t). This effect is further illus- 
trated in Fig. \7\ where the asymptote position is shown 
as a function of t together with yu (i) • Both curves are 
tangent at 74.91 /usee. Alternatively, the interpretation 
of this point can be understood from the inspection of 
Fig. Elb). In the (E , E)-p\a,ne this point is the intersec- 
tion between the LL and the E =4.906-eV line. On any 
constant-i curve the maximum value of y occurs in its in- 
tersection with the LL (point A), so in general this value 
of y is greater than the one at its intersection with the 
constant line E =4.906 eV (point B). But at the precise 
time where both lines intersect points A and B coincide. 
This is manifested in Fig. [7| where the y values of the 
asymptote position are always greater than those of the 
constant line E =4.906 eV except at the tangency point. 
Besides the mentioned asymptotic behavior, there ex- 
ist two other important features not described in the CA 
framework, viz. the main peak width and height de- 
pends on f, at variance with the left hand side member 
of Eq. (|25|l . Similarly, those behaviors are observed in 
deuterium (not shown in this paper). 



C. Study of the mean kinetic energy 

An immediate consequence of the present analysis can 
be found in the study of the mean kinetic energies, which 
is the main aim of the DINS technique. We will focus 
our attention to a sample constituted by an ideal gas of 
3 He with a mean kinetic energy of 30 K (ksT =1.7234 
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FIG. 6: (Color online) (a) Full lines: f(y,i) at 9 = 70° for 
the indicated times of flight. Vertical dotted lines: position of 
the asymptotes of f(y,t). Dotted lines: Resolution function 
times the amplitude factor £(£), centered at the y~M(t) values. 
See text for details, (b) Detail of the upper frame. Full thick 
line: f(y,t) for the special case of t=74.91 /xs, for which the 
asymptote position agrees with the value of y where the reso- 
lution function times £(t) (dotted thick line) is centered. Full 
thin lines: f(y,t) where the values of t are indicated in psec. 



meV). Such system was extensively examined in the 
Iiterature22i22i2i. The NCP corresponding to this system 
can be obtained in the exact formalism with Eq. i|21|) 
and is shown in Fig. EJa) (circles), which corresponds 
to a Gaussian distribution J(y) with a full-width-at-half- 
maximum (FWHM) of 1.3151 A -1 . The geometry, the 
filter and the efficiency are the same described at the be- 
ginning of this section. The function J(y) is shown Fig. 
lib) (dotted line). 

We now compare the exact formalism with the results 
produced in the CA framework. The resolution function 
can be obtained from the limit T — > OK of the NCP 
(shown in full line in Fig. |S{a)), and it is calculated in 
Eq. l(T2l . thus obtaining the Ri(y) shown in Fig. |8fb) 
(full line). The total cross section for gold employed in 
its calculation was taken from Ref. |27J. Following the 
procedure already described^ ., Eq. (|10f) was employed 
to perform a least-squares fit of a Gaussian J(y) on the 
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E= 4.906 eV 
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FIG. 7: (Color online) Asymptote position of f(y, t) as a func- 
tion of the time of flight t, for the case of hydrogen at a scat- 
tering angle of 9 = 70°, and flight length of Ln = 1105.5 cm 
and L\ = 69 cm. 



CA 



Mean Kinetic Energy (meV) 



Ri 
R> 
Rs 
Ri 



1.694 
1.374 
1.353 
1.519 



Exact value 



1.7234 



TABLE I: Mean kinetic energy of 3 He at 30 K obtained in 
the CA framework according to different resolution functions 
described in the text. The exact value is includes for compar- 
ison. 



exact NCP, letting the area and the width of the Gaussian 
as free parameters. Despite being R\ the only resolution 
function in the CA framework compatible with the limit 
at T ^0 K, we tested three other different resolution 
functions based on the customary use 

• i?2 (y) i obtained by fitting a Lorentzian function in 
the energy scal o 22 i 32 to the main resonance peak 
at 4.906 cV of the gold total cross sectional, and 
employing Eq (|12fl . 

• i?3(y), a Lorentzian function in the y variabl o 5 ' 6 ' 8 
fitted to R 2 (y)- The resulting FWHM is 3.292 A -1 . 

• Ri(y), a Lorentzian function in the y variable 568 
fitted to Ri(y). The resulting FWHM is 3.170 A" 1 . 

In table[I]we show the resulting effective temperatures, 
compared with the exact value, used as input to generate 
J(y). It is worth to remark that the different choices of 
R(y) produces large variations in the effective temper- 
ature, even when the FWHM has no such large varia- 
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FIG. 8: (Color online) (a) Neutron Compton profiles in the 
time of flight scale of 3 He at a scattering angle of 70 ° , at OK 
(line) and 30 K (circles), (b) The dotted line indicates the 
Gaussian momentum distribution J(y) of 3 He with a mean 
kinetic energy of 30 K. In full line, the resolution function 
calculated according to Eg. 11211 and c(€)t=ok shown in frame 
(a). Circles indicate the NCP calculated at 30 K expressed as 
a function of the y variable in the CA framework. 



tions. For example, the change in FWHM from R3 to 
i?4 is 3.7 % and the resulting mean kinetic energies differ 
in 12.3 %. The reason for such a large change, was re- 
ported in Ref. |2(J at low temperatures the shape of the 
observed NCP is determined mainly by the filter, so any 
imperfection in the description of the filter effects will 
affect significantly the obtained mean kinetic energy. 



V. CONCLUSIONS 

In this paper we developed the exact formalism that 
describes the NCPs obtained from DINS experiments in 
terms of an integral equation in the variable y that in- 
volves a kernel dependent on the experimental setup pa- 
rameters. For its calculation, a precise knowledge of the 
detectors efficiency, the filter total cross section and the 
incident neutron spectrum in a wide range of energies is 
required. As a result, we obtained a Volterra equation of 



the first kind that allows to obtain the desired momen- 
tum distribution J(y) by means of numerical methods. 

A significant part of this paper was devoted to the 
study of the variable y, that in the exact formalism here 
presented has a definition that differs significantly from 
that employed in the CA framework. The physical sig- 
nificance of this variable is the projection of the initial 
impulse of the target nuclei over the direction of the mo- 
mentum transfer vector q. In DINS data processing pro- 
cedures the passage from t to y is customarily calculated 
taking a fixed final energy, and thus it has a single well- 
defined value for each time channel. In this paper it was 
calculated exactly. The importance of knowing the exact 
mapping of the variable y in the (Eq, E^-plane can be 
referred to our recent study on the neutron final energy 
distributions in DINS experiments 26 . In that work we 
found that such distributions are far more complex than 
considering a single final energy, and depends on the time 
channel and the dynamics of the scattering species. We 
showed the constant-?/ curves on the (Eq , E)-plane for 
hydrogen and deuterium at a particular scattering angle 
and ratio of flight paths, which served to illustrate its be- 
havior, and developed its general expressions. A notice- 
able particularity of the y variable in the exact treatment 
that we present in this paper, is that it can have local 
extrema as a function of Eq for a given time of flight 
as shown in Fig. [21 which define straight lines in the 
(Eq, E)-pla,ne (LL). This causes singularities in the Ja- 
cobian dE /dyM which affect the kernel /m(2/m,£)- For 
example in the analyzed case of hydrogen, values of y 
greater that that those defined by the LL are not physi- 
cally allowed. 

In past publications^*^ we showed that CA can not 
be deduced from the exact expression. In this paper we 
again showed that the CA formalism is incompatible with 
the exact one, and that the basis of the difference between 
both formalisms resides in the integration kernels. Thus, 
the resolution function must be replaced by a kernel that 
depends on t. In Fig. I^Jwe show explicitly the differences 
between the exact kernel and the commonly accepted res- 
olution function. Besides the clearly observed differences 
in shape, it is notable the appearance of the above men- 
tioned singularities. In the particular case analyzed in 
this paper these discrepancies are particularly noticeable 
at 74.91 /xsec, where the singularity occurs exactly at the 
value of y where the maximum of the resolution function 
is centered, i.e. at y = yM(t)- 

So far, we have focussed our discussion on the differ- 
ences between the kernel /m (j/a/ , t) and the resolution 
function, without making any reference to the momen- 
tum distribution of the sample Jm^m)- As a first exam 
on the interrelation between fm(yM,t) and Jm(vm) it is 
helpful to check to what extent the singularity can affect 
significant parts of the integrand (|21|l . To this end we 
will compare the ratio E/Eq where the singularities ap- 
pear (defined by the slope of the LL) with the same ratio 
for the scattering on a sample defined by an ideal gas at 



10 



T=0K, i.e 



20 



E _ I cosO + ^/cos 2 6+{A-l)(A + l) 



E 



A+l 



(26) 



Eq. I|26|) defines also the most probable energy ratio in 
the case of a sample at finite temperature. In Fig. 0]we 
show in open symbols the ratio E/Eq of the LL for H, 
D and O, and in full symbols that corresponding to Eq. 
(|26[) . as a function of the scattering angle. We observe 
that for this particular value of I the singularities do not 
appear significantly near the main peak. However, due to 
the thermal motion (described by the shape of J(y)) the 
singularities affect the description of the tails of NCPs. 

The present formalism was formulated regardless of the 
function J(y). Although in this paper the case of an 
isotropic sample was examined, all the considerations we 
developed are applicable also for the case of anisotropic 
samples. It is also worth to notice that although in most 
of the significant cases concerning the study of Con- 
densed Matter the function J{y) is closely represented 
by a Gaussian shaped, this function could represent 
an anisotropic momentum distribution, not centered in 
y = like in the case of particles flowing in a prefer- 
ential direction. With regard to the analysis of experi- 
ments where the sample can be represented by the usual 
Gaussian distributions, it must be emphasized that the 
effects of the singularities will normally be smoothed out 
in the NCP, i.e. the experimentally accessible magni- 
tude. However, the differences between the exact for- 
malism and the CA in the description of the NCP are 
clearly manifested, as accounted in detail in Ref. l20t so 
we will briefly summarize them here. A comparison of 
the results produced by Eq. (|10fl in the calculation of 
the NCPs, with the exact formalism represented either 
by Eqs. Q) or (|2U[I (both formulations are equivalent, 
since they differ only in the integration variable), reveal 
that the CA is defective in the description of 

• the peak areas, 

• the peak centroids, 

• the peak widths, 

being the case of light scatterers at low temperatures 
the most sensitive ones to these imperfections. In this 
case the widths of the momentum distribution function 
J(y), and the resolution function R(y) becomes compara- 
ble, and thus the consequences of the differences between 
R(y) and the exact kernel is more clearly manifested on 
the NCPs. Thus, the obtainment of peak areas and effec- 
tive temperatures is significantly affected because of the 
use of the CA. The analysis shown in Sect. IIVCI clearly 
reflects this assertion, and the inadequacy to employ res- 
olution functions (either theoretically or experimentally 
defined) that due to the approximate nature of the CA 
has an unavoidable degree of ambiguity in its definition. 



We wish to conclude this paper giving a schematic out- 
line of the procedure we recommend to use to the exper- 
imentalists in the DINS field. For the measurement of 
an impulse distribution function it is necessary to char- 
acterize different parameters of the experimental config- 
uration. 

1. The incident neutron energy spectrum can be de- 
termined by a measurement with a detector of well- 
characterized efficiency placed on the direct beam. 

2. The efficiency of the detectors bank can be charac- 
terized using a heavy target (like Pb or Bi) at the 
sample position (with the filter out of the scattered 
neutron path). 

3. The effective filter thickness can be measured with 
a lead target, employing the method 'filter in - filter 
out', and fitting it as a parameter in Eq. Q). Al- 
ternatively, a transmission experiment can be per- 
formed on the direct beam. In the process, the full 
total cross section of the filter is needed. A wide 
variety of cross sections can be found v. gr. in Ref. 

E3 

4. Experimental data must be corrected by multi- 
ple scattering and attenuation effects, using Monte 
Carlo procedures like those described in Refs. 19 
and HI 

5. From the steps 1 to 3, which aim to know the kernel 
f{y,t), and step 4 which corrects for finite sample 
effects, we obtain the necessary elements to state 
Eq. (|21|) . that can be solved either by numerical 
means or by fitting parameters if we have a pre- 
vious knowledge of the function J(y). In this last 
instance the fitting process can be made also di- 
rectly on the t scale employing Eq. Q without 
changing the integration variable from Eq to y. 

Finally, let us mention that the present formalism is 
also applicable to the method of double differences re- 
cently presented^, which consists in employing two fil- 
ters of different thicknesses. The use of the formalism 
presented in our paper will be most beneficial for the 
users community of these DINS techniques. 



ACKNOWLEDGEMENTS 

This work was supported by ANPCyT (Argentina) un- 
der Project PICT No. 03-4122. We also thank Fundacion 
Antorchas (Argentina) , and CLAF (Centro Latinoameri- 
cano de Fisica) for financial support. 



APPENDIX A: DERIVATION OF dy/dE 



From the definition of the y of Eq. J7J) we have 



dy 
dE 



dy dq 



dq dE 



dy duo 



duj dE 
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Mw k\ dq 
2 ) dE^ 



<T 



M duj 



t q dEo 



(Al) 



So, we need to calculate the derivatives of q and uj with 
respect to Eq at constant t. Form Eq. Q) we have 



du) 

~d~E a 



l_L_dE_ 

h V dE 



From the kinematic condition (0) we have 



dE 
dE^ 



Lo(E_ 



3/2 



Replacing Eq. (|A3f) in Eq. (|A2I) we have 



dui 
dEo~ 



1 1+ W£ 
h\ ^ Li \e 



3/2N 



(A2) 



(A3) 



(A4) 



APPENDIX B: DERIVATION OF E(E ,y) 

From the Eq. (|13fl for y (which is in terms of E and 
Eq), we want to solve it for E, considering y and Eq as 
fixed input values. Defining the variables 



x 
a 

f 



E 
M 2 

2my 

m 

17 



(Bl) 



after some algebra Eq. i|13|) can De transformed to a 
fourth-order polynomial in x, with real coefficients 



^b lX l = 



(B2) 



On the other hand, deriving the definition of q presented 
in Eq. (J5J we have 



Then 

dq 
~dEo 



t h 2 q 



dq 

dEv~ 



dE n 



L fE 



L\ \Eq 



dq dE 



dEdEn 




(A5) 



cos 8 



(A6) 



where we have employed Eq. (|A3|l . Finally, replacing Eq. 
({A"4"|) and Eq. (|X5|) in Eq. I|A1|I we obtain Eq. (fT5jl . 



where 



b Q = \a (1 - nY E Q - lj E 

&i = 2 [1 + 2a (/i - ^i 2 ) £ ] ^f~E~ cos9 

b 2 = -l-2a[l- M 2 -2^ 2 cos 2 6»]S (B3) 

fe,3 = — 4a/^ (1 + [i) \/Eo cos 

64 = Q'(1 + m) ■ 



Therefore, the real positive roots of Eq. (|B2|) define the 
values of E(Eo,y), which define the y-constant curves 
represented in Figs. and 03 
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